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This  paper  derives  a  feedback  control  f(t),  Bf{t)l  <  r,  r  >  0,  which 

E 

forces  the  infinite  dimensional  control  system  =  Au  +  Bf ,  u(0)  =  u0  c  H 
to  have  the  asymptotic  behavior  u(t)  +  0  as  t-*-®  in  H.  Here  A  is  the 
infinitesimal  generator  of  a  CQ  semigroup  of  contractions  eAt  on  a  real 
Hilbert  space  H  and  B  is  a  bounded  linear  operator  mapping  a  Hilbert  space 
of  controls  E  into  H.  An  application  to  the  boundary  feedback  control  of  a 
vibrating  beam  is  provided  in  detail  and  an  application  to  the  stabilization 
of  the  NASA  Spacecraft  Control  Laboratory  is  sketched. 
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1.  Introduction 


/ 


^  In  this  paper '.we  consider^  the  feedback  stablization  of  a  linear  control 

•  A 

system  in  an  infinite  dimensional  state  space.  However  unlike  the  standard 

feedback  control  problem  where  the  goal  is  to  find  a  linear  feedback  control 

Clf 

law,  we  restrict  ourselves  to  the  case  where  the  controls  f(t)  satisfy  fche'if 
priori  constraint,  lf(t)B_  <  r,  r  >  0.  (Here  the  controls  f(t)  lie  in  a 

Ct 


and  |»|_  denotes  the  norm  of  E. )  This  contraint 
E 


Hilbert  space 

necessitates  a  choice  of  a  nonlinear  feedback  law  which  drives  our  state  u(t) 


to  zero  as  t  •. 

-»7<t  ♦••Sy  ^  -■ 

We  wi-Hr  derivei.such  a  nonlinear  feedback  law  based  on  "energy"  stability 

A 

methods.  The  analysis  of  the  asymptotic  behavior  of  the  state  u(t)  is  based 
on  the  theories  of  nonlinear  evolution  equations  and  contraction  semigroups. 
While  an  earlier  papers  {4-h  treated  a  related  problem  of  sub-optimal  control 
the  results  given  here  on  feedback  stabilization  are  new.  A  related  optimal 
control  problem  was  considered  by  Barbu.  [21  ^  ^ 

A  strong  motivation  for  this  paper  has  been  the  work  of  Hubbard  and  his 
co-workers  [3],  [4]  and  Balakrishnan  [5],  [6].  In  (3],  [4]  both  laboratory 
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experiments  and  computer  simulations  for  the  boundary  feedback  control  of  a 
vibrating  beam  were  described.  In  particular  the  problem  of  [3],  [4] 
corresponds  to  a  special  case  of  theory  given  here.  Hence  it  is  natural  to 
apply  the  abstract  results  of  this  paper  to  the  concrete  example  of  Hubbard  et 
al.  This  is  done  in  the  last  section  where  it  is  shown  that  the  feedback 
control  given  here  yields  a  stablizing  feedback.  Similarly  the  papers  of 
Balakrishnan  [5],  [6]  provide  a  mathematical  framework  for  the  stabilization 
problem  of  the  NASA  Spacecraft  Control  Laboratory  Experiment  (SCOLE).  The 
last  section  will  also  sketch  the  application  of  the  theory  given  here  to  that 
problem. 

The  paper  is  divided  into  seven  sections  after  this  one.  Section  2 
provides  a  statement  of  the  abstract  control  problem  and  a  hint  at  the  method 
of  resolution.  Section  3  gives  some  brief  preliminary  results  on  the  theory 
of  dynamical  systems  and  nonlinear  semigroups.  Section  4  uses  the  ideas  of 
Section  3  to  exposite  a  theorem  of  Ball  and  Slemrod  on  the  asymptotic  behavior 
of  a  class  of  nonlinear  evolution  equations.  Section  5  applies  this  theorem 
to  yield  one  resolution  of  the  feedback  stabilization  problem  (Theorem  5.1). 
Section  6  gives  a  survey  of  the  results  on  asynptotic  behavior  of  nonlinear 
contraction  semigroups  and  a  useful  theorem  of  Dafermos  and  Slemrod  is 
presented.  Section  7  applies  Dafermos  and  Slemrod's  result  to  the  feedback 
stabilization  problem  in  the  case  A  has  compact  resolvent  and  E  =  R 
(Theorem  7.1).  Section  8  uses  the  earlier  mentioned  problems  of  Hubbard  et  al 
and  Balakrishnan  as  illustrative  examples. 

We  note  that  a  good  reference  for  the  ideas  on  nonlinear  semigroups  and 
asymptotic  behavior  is  the  monograph  of  A.  Haraux  [7].  Many  of  the 
propositions  used  here  may  be  found  there.  In  addition  numerous  examples 
illustrating  the  nonlinear  semigroup  theory  are  contained  there  as  well. 


» 


2.  The  control  problem 

Let  H  be  a  real  Hilbert  space  with  inner  product  <*,•>  and  norm  |*|  . 
and  let  E  be  a  second  real  Hilbert  space  with  inner  product  (•»•)  and  norm 
|«IE.  Also  let  A  be  the  infinitesimal  generator  of  a  linear  CQ  semigroup 
of  contractions  on  H  denoted  by  eA  .  (In  particular  we  know  |e  u^l  < 
luQl  for  all  u0  €  H  and  <A$,$>  <  0  for  all  4>  «  D(A).)  Finally  let  B 
be  a  bounded  linear  operator  from  E  to  H. 

We  consider  the  abstract  control  system 

^7-  -  Au  +  Bf  ,  (2.1) 

dt 

u(0)  -  u0  €  H  .  (2.2) 

Our  goal  is  to  find  a  feedback  control 

f(t)  -  K(u(t) )  (2.3) 

satisfying  the  constraint 

lflE  <  r,  r  >  0  ,  (2.4) 

which  will  yield  u  -  0  globally  asymptotically  stable  in  some  sense. 

To  do  this  first  formally  coopute  the  time  rate  of  change  of  the 
"energy": 


lu(t)  I*  “  <Au,u>  +  <Bf,u> 
2  at  n 


(2.5) 


<  (f,B  u) 

Here  we  have  used  the  fact  <Au,u>  <  0.  In  order  to  force  energy  decay  yet 
satisfy  the  constraint  (2.3)  we  use  a  saturating  control  law  as  suggested  in 
the  work  of  Gutman  [8]  for  finite  dimensional  systems,  i.e.  we  set 


« 

K(u)  -  -  if  IB*ulE  >  r 

IB  UIE 


■  -B  u 


if  IB*ulE  <  r 


(2.6) 


-3 


Notice  K(u)  is  continuous  as  a  function  of  u  and  that  the  desired  energy 
dissipation  is  obtained: 


4  *u<t)|2  <  -rlB*uI  if  |B*ul  >  r 

<  at  E  E 


*  2  * 

C  -IB  ulE  if  IB  ulE  <  r  . 


(2.7) 


In  the  next  section  we  will  analyze  (2.1)  with  feedback  f(t)  *  K(u(t)) 


and  give  sufficient  conditions  for  its  successful  implementation. 


Preliminary  results  on  nonlinear  semigroups 


Definitions.  Let  H  be  a  real  Hilbert  space.  A  (generally  nonlinear) 
semigroup  T(t)  on  H  is  a  family  of  continuous  maps  T(t)  :  H  +  H,  t  e  R+, 
satisfying  (i)  T(0)  *  identity,  (ii)  T(t+s)  *  T(t)T(s),  for  all  t,  s  e  R+ . 
If  in  addition  IT(t)$  -  T(t)i|il  <  l^-ijil  for  all  $,  c  H,  t  >  0,  T(t)  is 


called  a  contraction  semigroup. 

For  4  e  H  define  the  positive  orbit  through  $  by  0+($)  *  U  T(t)<ji. 

tcR 

The  ui-limit  set  of  <p  is  the  (possibly  enpty)  set  u»(^)  *  {4»  e  H;  there 


exists  a  sequence  tn  -*■  <»  as  n  ♦  «•  such  that  T(tn)$  +•  ifi  as  n  >  *} .  The 
weak  ai-limit  set  of  $  is  the  (possibly  en^ty)  set  given  by  e  H; 

there  exists  a  sequence  tR  •  as  n  -*•  <*>  such  that  T(tn)$-k  i|>  as  n  *  «•} . 
(Here  the  symbol  denotes  weak  convergence  in  H.) 


As  subset  C  of  H  is  said  to  be  positively  invariant  if  T(t)C  C  C 
for  all  t  e  R+,  and  invariant  if  T(t)C  =  C  for  all  t  e  R+. 


Theorem  3.1.  (i)  If  0($)  is  precompact  .£hen  o>((f>)  is  a  nonenpty,  invariant 

set  in  H.  (ii)  If  each  T(t)  is  sequentially  weakly  continuous  on  H  (i.e. 
T(t)$n  -*  T(t)$  is  ♦),  then  ()  +  ($)  bounded  implies  is  a 

nonempty,  invariant  set  in  H. 


ram»uiTiwmwwwmjpiKHmBt:<»i«TmBMiinMCTiwM.miihnuw»wimwHinirmnnrHwu»»ruirvirB 


Proof .  (i)  The  proof  is  a  direct  consequence  of  Prop.  2.2  in  Dafermos  [9]. 

(ii)  Since  (J*"( <^)  belongs  to  a  sequentially  weakly  compact  set  in  H,  ojw($) 
is  non-en$>ty.  Furthermore,  since  H  is  separable  this  weakly  compact  set  may 
be  regarded  as  a  compact  set  in  a  metric  space  induced  by  the  weak  topology 
(see  Dunford  and  Schwartz  [10]).  The  result  again  follows  from  Prop.  2.2  in 
Dafermos  [9]. 

Hidden  in  Theorem  3.1  is  the  essence  of  this  paper.  Namely  that  in  the 
study  of  the  feedback  control  system  described  in  Section  2  we  may  need  to 
use,  under  different  circumstances,  either  part  (i)  or  part  (ii)  of  Theorem 
3.1.  Roughly  the  idea  is  that  in  the  study  of  nonlinear  semigroups  of 
"parabolic"  type  and  nonlinear  contraction  semigroups  of  "hyperbolic"  type 
sufficient  conditions  have  been  given  for  0+($)  to  be  precompact  and  hence 
<ii(<|i)  to  be  nonempty  (see  Henry  [11],  Pazy  [12],  and  Dafermos  and  Slemrod 
[13] J.  On  the  other  hand  other  applications  may  yield  only  the  information 
that  0+(<fr)  is  bounded  and  hence  the  main  tool  in  studying  the  asymptotic 
behavior  of  the  feedback  system  will  be  the  weak  ui-limit  set. 


4.  Semi linear  evolution  equations 

We  recall  some  standard  results  on  nonlinear  evolution  equations. 
Consider  the  initial  value  problem 


du 

dt 


Au(t)  +  F(u(t) ,t)  , 


u(t0) 


I 


(4.1) 


where  A  is  the  infinitesimal  generator  of  a  linear  Cg  semigroup  eAt  on  a 
real  Hilbert  space  H  with  inner  product  <•,•>  and  norm  II,  F  :  H  x  R  ♦  H 
is  a  given  function  and  ug  «  H  is  a  given  initial  datum. 
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Definition.  Let  tj  >  tg.  A  function  u  e  C([tQ,t^]?H)  is  a  weak  solution 
of  (4.1)  on  [tg,t.,]  if  u(tQ)  -  uQ,  F(u( •  )  *  •  )  t  L 1  ( t o / 1 1  j H )  and  if  for  each 
w  «  D(A*)  the  function  <u(t),w>  is  absolutely  continuous  on  [tg,t.j]  and 
satisfies 

^—<u(t),w>  *  <u(t),A*w>  +  <F(u(t) ,t) ,w> 
at 

for  almost  all  t  e  [tg,tj]. 

Theorem  4.1.  (cf.  Balakrishnan  [14],  Ball  [15]).  Let  t^  >  tg.  A  function 
u  :  [t0,t.|]  ♦  H  is  a  weak  solution  of  (4.1)  if  and  only  if  F(u(t),«)  e 
L^  (tQ,t<j  ;H)  and  u  satisfies  the  variation  of  constants  formula 

u ( t )  =  e  0  uq  +  /*  eA^t“s^F(u(s) ,s)ds 
for  all  t  c  [  t  g ,  t  ^  ] • 

The  next  result  characterizes  the  asymptotic  behavior  of  solutions  to 
(4.1)  in  an  important  special  case.  Also  we  assume  system  (3.1)  is 
autonomous,  i.e.  F(t,u)  *  F(u),  tg  *  0. 

Theorem  4.2.  Let  A  generate  a  linear  Cg  semigroup  eAt  of  contractions. 
Let  F  :  H  +  H  satisfy 

(i)  F  is  locally  Lipschitz 
(ii)  4»n  -*>  ill  ««>  F ( )  ♦  F( 4»)  , 

(iii)  <F(\Ji)  ,<j»>  <  0  for  all  i|>  c  H. 

Then  (4.1)  possesses  a  unique  weak  solution  u(t;ug)  on  R+  for  each  Ug  e  H 
Furthermore  T(t)ug  *  u(t,ug)  defines  a  semigroup  on  H,  ^(ug)  is  a 
nonempty  invariant  set  for  each  Ug  e  H,  and  for  each  c  u^Cug) 

<T(t)<i,  F(T(t)iJ)>  -  0  for  all  t  c  K+  • 

If  in  addition,  the  only  solution  to  the  above  equation  is  <j;  *  0,  then 
u(t;ug)  -*0  as  t  ♦ 

-6- 
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Proof,  foie  proof  is  given  in  the  paper  of  Ball  and  Slemrod  [16].  A  central 

idea  of  the  proof  is  the  dissipative  mechanism  (iii)  yields  0+(uq)  bounded 

and  hence  by  Theorem  3.1  (ii)  ui^Uq)  is  a  nonempty,  invariant  set  in  H. 

2 

The  Liapunov  functional  flu ( t )  H  is  then  used  to  identify  co^Uq)  as  noted 
in  the  theorem.  v 


5.  Application  of  Theorem  4.2  to  the  stabilization  problem 

In  this  section  we  will  discuss  the  asymptotic  behavior  of  the  feedback 
system  (2.1),  (2.2),  (2.5),  (2.6),  i.e.  we  will  study  the  nonlinear 
evolutionary  system 


du 

dt 


Au  +  G(u) 


(5.1) 


u(0)  =  u0 


(5.2) 


with 


* 

rBB  n  * 

G(u)  =  -  — if  IB  u)  >  r 


IB*ul 


E 


*  * 

-BB  u  if  IB  u|  <  r 


(5.3) 


where  A,  B 


are  as  given  in  Section  2. 


Theorem  5.1.  For  each  Uq  e  H  there  is  a  unique  weak  solution  u(t;u0)  = 
T(t )Up  of  (5.1),  (5.2)  defined  for  all  t  e  R+  with  {0}  a  stable 
equilibrium.  If  in  addition  B  is  conpact  and  the  only  solution  of  the 


equation 

B*eAt  ij)  =  0  for  all  t  e  R+ 

is  \|)  s  0,  then  u(t;ug)  0  as  t  -*■  <»  for  all  u0  e  H* 


(5.4) 


-7- 


Proof.  First  we  note  G(u)  is  globally  Lipschitz  continuous.  For  if  u,, 
u2  €  H  and  we  set  y,  -  B*uv  y2  -  B*u2  then  ly^  <  r  and  |y  »  >  r 


implies 


■ac«1>-ccu2>i  <  II | b i | |  iYl  -  Tp-f-ij 


<  r  ,yi,y2flE  '  ry2"E 


<  ^  '*i<>*2'e  -  ly,!,)  +  7,17,1.  -  ry2.E 

<  iiMli  («yilE  fly2-yi»E  +  ly,(ly,lB-r)  +  rCy^)^ 

<  II lBl I  I ( |72-y1 «E  +  I ■y1 »  “  r|  +  ly2~y1 g ) 

<  3I I  lBl  1 1  ly^y^  • 

Here  we  have  used  |  |  |  •  |  |  |  to  denote  the  operator  norm  on  the  space  of 
bounded  linear  operators  E  H.  On  the  other  hand  when  both  y y2  have 
norm  greater  than  or  equal  to  r 

|G(u  )-G(u  )|  <  |  |  | B f  |  | r  n- J— - — g 

2  iy,*B  »y2iB  * 

‘  iy,i*iyJE  ,yi,y2,E  -  We'e 

*  r  'We "  +  -  We'e 

« ilisill  (2lyi,s  ,VVe) 


<  2 1 1 1 B 1 1 1  ,y2-y1,E  • 

Of  course  if  ly,!  <  r.  iy^l  <  r  then  IGlu,)  -  GIUjH  <  |||b|||  ,^-y^  so 

in  general  W(U) )  -  GlUj) I  <  3| | | b| | | 2|u,-u2l  for  .11  u,,  u2cH. 

Next  note  that  the  compactness  of  B  implies  G(iJ,  )  +  G(0»)  if  a, 

n  *n 

Also  <G(if>)  ,$>  -  -rlB%lE  if  lB%l>r  and  <G(*),*>  -  -IB%g*  if  gB%|  <  r. 


-8- 


So  Theorem  4.2  applies  and  tells  us  that  for  each  uQ  c  H  (^(ug)  is  a 
nonempty  invariant  set  in  H  and  for  each  c  uw^u0^ 

<T(t)i|»,G{T(t)4»)>  =  0  for  all  t  e  R+  .  (5.5) 

But  <T(t)i|),G(T(t)i|)>  =  0  for  all  t  e  R+  implies  B*T(t)4>  *  0  for  all 
t  e  R+  which  in  turn  tells  us  G(T(t)<}>)  =  0  for  all  t  e  R+.  But  then  the 
variation  of  constants  formula  of  Theorem  4.1  shows  for  i(i  e  m^Ug)  that 
T(t)iJ>  =  eAt^-  So  (5.5)  in  fact  implies  (5.4)  and  hence  tiiw(ug)  =  {0}.  The 
fact  that  {0}  is  stable  trivially  follows  from  the  estimate 
IT(t)uQ I  <  JuQ(. 


Lipschitzian,  and  dissipative  i.e.  <Q(u),u>  <  0  for  all  u  e  H.  Then  the 
conclusion  of  Theorem  5.1  still  holds  where  f(t)  =  K(u(t)). 


Proof .  If  we  insert  f(t)  =  K(u(t))  into  (5.6)  our  feedback  system  is 

*  Au  +  Q(u)  +  G(u) 
dt 

But  now  that  argument  given  in  the  proof  of  Theorem  5. 1  applies  verbatum 
with  G  replaced  by  Q  +  G. 


6.  Asymptotic  behavior  of  nonlinear  contraction  semigroups. 


Let  T(t)  be  a  nonlinear  semigroup  of  contractions  on  a  real  Hilbert 


space  H.  We  denote  D (A)  be  the  set  of  those  <(>  e  H  for  which 


lim 
h>0  + 


T(h)j)-il) 

h 


exists  and  define 


-A*  =  lim 

h+0+  h 

It  is  well  known  from  the  theory  of  nonlinear  contraction  semigroups 
(e.g.  [17],  [7])  that  associated  with  a  nonlinear  semigroup  of  contractions 
there  is  a  unique  (possible  multi-valued)  operator  -A,  that  A  is  maximal 
monotone,  D(A)  is  dense  in  H,  range  (XA+I)  =  H  for  any  \  >  0,  and 
(\A+I)_1  is  a  continuous  single-valued  function. 

On  the  other  hand  given  a  maximal  monotone  operator  A  we  know  for  every 
Uq  c  D(A)  there  exists  one  and  only  one  function  u(t)  :  (0,®)  >  H  such  that 


u ( t )  e  D(A)  , 

du{  t ) 


dt 


c  L  [  (0 ,») ?  H]  with 


1§I  <  IA°V 


+  Au(t)  »  0  on  (0,®)  , 


du(t ) 
dt 

U(0)  =  Uq  • 

Here  A°,  the  miniirial  section  of  A •  is  the  function  which  assigns  to  each 
$  €  D(A)  that  element  of  A$  which  has  least  norm.  Furthermore  u  is  right 
-  differentiable  at  any  t  e  [0,®)  and 


for  all  t  e  [0,®). 

If  u  and  v  are  solutions  associated  to  initial  data  Uq  and  vq 

then 


lu(t)  -  v(t)lH  <  fluo”vollH  f°r  >  0  * 


-10- 
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We  note  by  T(t)  the  extension  by  continuity  of  the  map  Ug  e  D{A)  + 
u(t)  €  D(A)  to  D( A) .  If  D(A)  is  dense  in  H  then  T(t)  defines  the 
nonlinear  contraction  semigroup  on  H  "generated"  by  -A. 

The  asynptotic  behavior  of  nonlinear  contraction  semigroups  has  been 
characterized  by  the  following  theory  of  Dafermos  and  slemrod  [13]. 


Theorem  6.1.  Let  A  he  a  maximal  monotone  operator  on  a  Hilbert  space  H. 
Assume  0  e  range  (A)  and  (XA+I)”1  is  conpact  for  some  X  >  0.  Then  for 
any  uQ  e  D(A)  the  weak  solution  of  the  Cauchy  problem 


+  Au(t)  3 


0 


u(0)  =  uQ  , 

given  by  u(t)  *  T(t)u0  approaches  as  t  ®  a  conpact  subset  ft  of  a 

sphere  {y;  ly-al  =  r},  r  <  lu^-al,  a  e  A~*0.  Furthermore  ft  is  minimal, 

invariant,  and  equi-almost  periodic  under  the  semigroup  T(t)  generated  by  -A 

and  T  restricted  to  the  closed  convex  hull  of  ft  coft  is  an  affine  group  of 

isometries.  If  in  addition  ug  €  D(A)  then  the  set  ft  is  contained  in 

D(A),  and  A°ft  is  conpact  and  lies  on  a  sphere  centered  at  0.  Moreover 

coft  C  D(A)  and  the  restriction  of  A®  to  coft  is  affine. 

While  the  proof  is  contained  in  [13]  and  [7]  we  note  the  main  idea  again 

relates  back  to  Theorem  3.1.  Recall  that  we  noted  after  the  statement  cf 

Theorem  3.1  that  part  (i)  of  the  theorem  applied  to  certain  "hyperbolic" 

semigroups  where  we  could  actually  show  0+($)  is  precompact.  Theorem  6.1 

identifies  a  class  of  these  semigroups  as  ones  arising  from  nonlinear 

contraction  semigroups  with  compact  resolvent.  The  rest  of  the  theorem 

essentially  follows  by  identifying  the  w-limit  of  T(t)$  with  the  aid  of  the 

2 

Liapunov  functional  lu(t)|  . 


k 

I 

t 

i 


-i  i- 


vw-. 
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slication  of  Theorem  6.1  to  the  stabilization  problem. 


In  this  section  we  discuss  the  application  of  Theorem  6.1  in  the  special 


case  E  *  R.  In  this  case  B  is  a  fixed  element  in  H.  Furthermore  note  . 
that  for  ij;  e  D(A)  that  -A-G  is  monotone.  To  see  this  set  y  =  B%» 
z  *  B%  and  observe 


(i)  for  lyl  <  r,  Izl  <  r  , 

E  E 

A4>  +  G(  4>)  -  AiJ;  -  G  ( -4» )  >  <  -|y-z|2? 

(ii)  for  lyl  >  r,  Izl  >  r  , 

E  E 

<$-!)»,  a{>  +  g( )  -  a<ji  -  G(t|») >  <  -r(y-£)-(y^j-  -  yjfy)  <  o  » 

(iii)  for  lyl  <  r,  Izl  >  r, 

E  E 

<tj>-<j>,A<|»  +  G(^>)  -  a*  -  G{(j») >  <  -(y-£)  (y  ■|||“) 

<  -y2  +  r|y|  +  I y |  Icl  -  dd 
( |y|-r)  (|d-|y| )  <  -<|y|-r)2. 

Bence  we  see  that  the  operator  A  *  -A  -  G  defined  on  D(A)  is  monotone. 
Notice  also  that  the  above  argument  also  shows  -G  is  monotone.  Since  -G 
is  continuous  (see  Section  5)  and  -A  is  maximal  monotone  (by  the  Hille- 
Yosida -Phil lips  Theorem)  a  theorem  of  G.  F.  Webb  [18]  asserts  that  the  sum 
-A  -G  is  maximal  monotone. 


Now  we  are  prepared  to  prove  our  stabilization  theorem. 


Theorem  7.1.  For  each  uQ  e  H  there  exists  a  unique  weak  solution  of  (5.1), 
(5.2)  for  all  t  e  R+  with  {0}  a  stable  equilibrium  of  (5.1).  If  in 
addition  E  -  (XI-A)"1  is  compact  for  all  real  X  >  0,  and  the  only 

solution  of  the  equation 

B*eAti|>  ■  0  for  all  t  e  R+ 

is  $  ■  0,  then  u(t,ug)  ♦0  as  t  ♦  «•  for  all  Uq  €  H. 

Notice  here  that  Theorem  7.1  improves  on  Theorem  5.1  in  that  weak 
convergence  is  now  replaced  by  strong  convergence.  Of  course  the  price  paid 


is  that  we  assume  E  *  R  and  A  has  contact  resolvent. 

Proof  of  Theorem  7.1.  We  shall  apply  Theorem  6.1.  We  have  already  shown 
»  -A  -  G  is  maximal  monotone  and  also  trivially  0  «  range(A).  Now  let 
{1^}  be  a  bounded  sequence  in  H  with  e  D(A)  such  that 

u  h 

-»“» -  G<v  ♦  f  -  r  ,,-1) 

for  X  >  0.  Take  the  inner  produce  of  both  sides  of  (7.1)  with  un  to  see 

lu  |2  <  <u  ,h  >  <  (u  II  |h  B 
n  n  n  n  n 

and  hence  lu^l  <  ,hn**  Thus  {v^}  also  belongs  to  a  bounded  set  in  H. 
Rewrite  (7.1)  as 

1  h 

-Ta-vDu  -  ■  +  G(u  )  .  (7.2) 

v  X  n  X  n 

Hie  Lipschitz  continuity  of  G  shows  the  right  hand  side  of  (7.2)  lies  in  a 
bounded  set  of  H  and  the  compactness  of  (A-XI)-1  for  X  >  0  shows  {ufl} 
lies  in  a  compact  subset  of  H.  Thus  (X  +I)-1  is  compact  for  all  X  >  0. 

Now  let  Uq  e  D(A) .  Theorem  6.1  tells  us  u(t)  =  T(t)ufl,  the  strong 
solution  of  (5.1),  (5.2),  approaches  as  t  ■*>  »,  a  contact  subset  0  of  a 
sphere  {yj  lyl  =  r}  where  0  C  D(A).  Let  Vq  e  fi.  Since  £)  is  invariant 
we  must  have  HT(t)Vgl  =  r  for  all  t  e  R+,  and  differentiation  with  respect 
to  t  and  use  of  (5.1)  shows  B*T(t)Vg  =  0  for  t  e  R+ •  But  (7.4)  coupled 
with  (5.1)  and  the  invariance  of  (1  shows  T(t)vQ  =  ®  v0  for  vQ  c  ft. 

Hence  for  vQ  c  fl,  B*eAtvQ  ■  0  for  all  t  «  R+.  But  by  the  hypothesis  of  the 
theorem  Vg  »  0  and  Q  ■  {0}  for  Ug  e  D(A).  Since  D(A)  is  dense  in  H 
and  T(t)  is  a  contraction  the  triangle  inequality  readily  shows  £2  *  {0} 
for  all  Ug  c  H  as  well.  Since  the  existence,  uniqueness,  and  stability  have 
already  been  prove  in  Theorem  5.1  the  proof  is  complete. 
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where  we  require  a(t)  =  — 

at 


,  b(t) 


x  *  L 


x  =  L 


The  above  first  order  formulation  motivates  us  to  formally  define  the 


linear  operator  A 


dx  ] x  =  L 


2  x  =  L 


so  that  when 


,  B  - 


(5.3)  has  the  form 


Au  +  Bf  • 


(8.4) 


To  be  precise  we  must  identify  the  domain  of  definition  of  A.  To  this 


extent  we  first  define  the  Hilbert  space  H: 


H  -  {(w,V,a,b)  €  H2  (0,1)  x  L2(0,1)  x  R  x  R  ; 


w  -  w'  »  0  at  x  -  0} 


endowed  with  the  inner  product 


<(w,v,a,b),  (w,v,a,b)> 


/q  (w''(x)  w''(x)  +  v(x)v(x))dx  +  aa  +  bb 


Now  we  take  the  domain  of  A  as 


D(A)  -  {(w,v,a,b)  €  H  (0,1)  x  H  (0,1)  x  R  x  R  ; 


av  * 

W  -  w.  -  o,  V  -  V  -  0  at  x  -  0,  V  =  a,  —  =  b  at  x  -  l} 


It  is  an  easy  computation  to  see  that  <Au,u>  is  dissipative  for 


u  c  D(A).  In  particular  this  shows  that  A  is  dissipative  i.e.  <Au,u>  <  0 


for  u  e  D( A)  or  equivalently  -A  is  monotone  <-Au,u>  >  0  for  u  e  D(A). 


We  also  note  that  A  is  an  infinitesimal  generator  of  a  linear  Cr 


semigroup  of  contractions  T(t)  on  H.  To  show  this  we  simply  apply  the 


iAimer-Phillips  Theorem  (18)  which  asserts  A  will  be  the  generator  of  such  a 


semigroup  if  and  only  if  A  is  dissipative,  densely  defined  on  H,  and 


satisfies  the  range  condition  R(-A  +  AnI)  *  H  for  some  An  >  0. 


We  have  already  shown  A  is  dissipative  and  it  is  a  simple  observation 


to  see  that  it  is  densely  defined.  To  check  the  range  condition  we  let 


(g,h,c,d)  be  a  generic  element  of  H.  Then  satisfying  the  range  condition  is 


equivalent  to  finding  u  -  (w,v,a,b)  e  D(A)  with 


-v  +  A0w  -  g 


(8.5) 


4 

d  w  .  ,  , 

— -  +  Aqv  -  h  , 

dx 


(8.6) 


.3 

d  w  |  .  , 

“a  ,  *  V  c  ' 

dx  1  x  ««  L 


(8.7) 


.2 

d  w  i 

~2  +*0*-*  ■ 
dx  ‘ x  *  L 


(8.8) 


From  the  first  two  of  these  equations  we  see  w  should  satisfy 


d^w  2  . 

—  +  Aqw  -  h  +  A0g 

dx 


(8.9) 


We  now  solve  the  ordinary  differential  equation  subject  to  the  boundary 


conditions  w(0)  ■  0,  w’(0)  *  0,  and 


,  S. 
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d  w  2 

—  +  XQw  =  C  +  XQg  , 
dx 


±Jt+x2±L=d  +  x  42 

„  2  A0  dx  0  dx 


(8.10) 


(8.11) 


at  x  *  L.  Thi3  can  be  done  by  explicitly  solving  (5.9)  with  four  constants 
of  integration  and  then  using  the  four  boundary  conditions  to  evaluate  the 
constants.  This  will  yield  w  e  H4(0,1)  with  w(0)  *  w'(0)  =  0,  which 
with  v  defined  by  (5.5)  and  a,  b  given  by  a  =  v(L),  b  *  (L)  solve 

(5.5)— ( 5.8) .  Straightforward  inspection  of  (5.5)-(5.8)  also  shows  u  = 
(w,v,a,b)  is  in  D(A) . 

This  analysis  works  for  any  Xq  real.  This  is  no  surprise  since  as  we 
shall  prove  the  spectrum  of  A  is  purely  imaginary,  discrete,  of  the  form 
X  «  ±  iy2  where  y  satisfies  the  transcendental  equation 


1  +  cosyL  coshyL  +  y(sinhyL  cosyL  -  coshyL  sinyL) 

3  4 

-  y  (cosyL  sinhyL  +  sinyL  coshyL)  +  y  (1  -  cosyL  coshyL)  *  0  . 


(8.12) 


Finally  we  trivially  note  that  B  is  a  bounded  linear  operator  on  the 
control  space  E  -  R  to  H.  Hence  we  have  rewritten  the  boundary  control 
system  (8.1),  (8.2)  in  the  form  (2.1),  (2.2).  Here  the  initial  data  is 


and  we  set  (a, a)  *  aa,  a, a  e  R. 

£ 

Also  since  the  imbedding  of  h4(0,1)  x  H2(0,1)  +  H2(0,1)  x  L2(0,1)  is 
compact  we  see  D(A)  is  conpactly  imbedded  in  H.  Hence  (A  -  Xgl)~^  is 
canpact  for  any  real  X«. 


I 


8 

* 

$1 

1 


s' 


I 
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1 
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So  far  we  have  shown  A  is  the  infinitesimal  generator  of  a  Cq  semi¬ 
group  of  contractions  T<t)  on  H,  (A-XqI)”1  :  H  -*•  H  is  conpact.  To  apply 
Theorem  7. 1  we  must  now  show  B*eAti|>  “0,  ip  e  H,  for  all  t  e  R+,  implies 
^  ■  0.  But  if  B*eAttJ>  *  0  then  <Bq,  eAt4>>  =  0  for  all  q  e  R.  But  by  the 
definition  of  B  if  we  write  u(t)  =  eA%  =  [w(t i  ,v(t)  ,a( t)  ,b(t)  ]  then 

A  A 

qb(t)  *  0  for  all  q  e  R  and  t  e  R ,  i.e.  b(t)  =  0  for  all  t  c  R  .  So 

*  + 

our  goal  now  is  to  show  b(t)  =  0  for  all  t  e  R  implies  i|i  *  0. 

To  do  this  we  shall  conpute  eAti(i  for  ip  e  D(A)  explicitly.  First  let 
us  extend  H  to  the  complex  Hilbert  space  H  m  H  ®  iH  where  H  has  inner 
product  <<  ,  >> 

<<Xj  +  iyi ,  x2  +  iy2>>  = 

<X1/X2>  +  <yi/V2>  “  1  t<xl'y2>  “  <X2'yi>3  * 


Then  A  is  also  the  infinitesimal  generator  a  Cq  semigroup  on 

1 1  At  At  At  v  —  ii 

H  :  e  xm«  He  x  +  ie  Im  x  ^or  X  €  n 

The  advantage  of  introducing  H  is  that  we  can  now  represent  eAtufl  as  an 

eigenfunction  expansion.  To  this  end  let  u  c  D(A)  so  that  Au  *  Xu.  Then 


if  we  write 


u  must  satisfy 


we  see  w  must  satisfy 
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£ 


at  x 


(8.15) 


w""(x)+Xw(x)=0,  0<x<L  ,  (8.13) 

w  m  w*  =0  at  x  =  0  ,  (8.14) 

-w*  * •  +  \2w  =  0 

w"  +  x2w'  =  0 

From  (8.13),  (8.14)  we  know  w(x)  is  of  the  form  w(x)  =  C(sin  px  - 
sinhpx)  +  D(cos  px  -  coshpx)  where  X  »  ±  ip  .  From  (8.15)  we  seen  the 
additional  relations 

C(sin  pL  +  sinh  pL)  +  D(cos  pL  +  cosh  pL)  = 

-p2{C(cos  pL  -  cosh  pL)  +  D (-sin  pL  -  sinh  pL)} 

and 


(8.16) 


(8.17) 


C(cos  pL  +  cosh  pL)  +  D(-sin  pL  +  sinh  pL)  = 
p{C(sin  pl>  -  sinh  pL)  +  D(cos  pL  -  cosh  pL)} 
must  be  satisfied. 

Elimination  of  C  and  D  from  (». 16) ,  (8.17)  yields  the  spectral 
formula  (8.12).  For  each  p  satisfying  (8.12)  the  associated  C  and  D  are 
determined  from  either  (8.14)  or  (8.17).  If  we  denote  the  positive  solutions 
of  (8.12)  by  {pn}  then  0<p1<p2<  ■*■+»  and  the  eigenvalues  are  XR  = 

2  „  «  1  T  "3  .oonr-i  plapnfunctions 


m 

2 

+lPjj,  n  *  1,2,...,  X_n  = 

s 

are 

fc 

r 

+ip2wna) 


+ip2wn(L) 


for  X, 


1,2. 


and 
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aaa 


fwhu 


kr,  TV*  ^  *•  "  K."  K  ”  ■JTK”  Kr  Kr  »L’  X."  *T  v"  *."  X' 


-ip^w 

n  n 


for  X_n  ,  n  =  1,2. . . 


-L<wn{L) 


-iUnwn(L> 


(Here  we  have  used  vn  =  Xn  wn* ) 


From  (8.13),  (8.14),  (8.15)  we  find 


w"(x)w"(x)dx  +  X2  w  (x)w  (x)dx  +  X2  [w  w  +w  w  ]  =0  .  (8.18) 

J  0  n  m  +n-'0n  m  +n  n  m  n  m  x=L 


Now  interchange  m  and  n  in  (8.18)  and  take  the  difference  of  the  two 


equations.  This  will  show  that 


,L 

|„  w  (x)w  (x)  +[ww+ww]  =0  , 

■'On  m  nm  nm  x=L 


fr  w"(x)w"(x)dx  -  0  , 

J  0  m  n 


if  m  f  n  . 


(8.19) 


If  we  set  m  *  n  (8.18)  we  see  as  well  that 


/  w"(x)  dx  =  y  {/  w  (x)dx  +  [w  +  w  ]  } 

J  o  in  n  0  n  n  «  v*t. * 


n  n  x*LJ 


(8.20) 


From  (8.19)  we  see  that  <<u  ,u  >>  =0  if  m  f  n  and  m  ^  n  and  from 

n  m  i  i  i  i 


(8.20)  we  see  <<un,u_n>>  =  0  as  we^1»  n  =  1,2,...  .  Hence  (un),  n  = 


±1,±2,...  of  eigenfunctions  forms  an  orthonomal  set  in  H.  It  is  easy  to 


show  it  is  complete  as  well.  Furthermore  we  assume  it  is  normalized  so  that 


<<un*un>>  *  1.  Hence  if  uQ  e  D(A)  C  H 


*  .  At  r  nt 

u(t)  “  e  u.  =  )  c  e  u 

0  *•  n  n 


(8.21) 


where  the  (cn)  are  complex  coefficients  satisfying 


cn  “  <<u(j,un>>,  n  *  ±1,  ±2,... 


(8.22) 


V.V.S  '■  .•> .%  „■»  ■*  .%  .>  .V  A  .V  >  .y*w 


Since  u  =  u  .X  =  X  ,  (where  the  overbar  denotes  conrolex 
n  -n  n  -n  ^ 


conjugation)  we  can  write  (8.21)  as 

*  "  Xnt  “  Xnt 

u(t)  =  Tce  u  +  Tc  e  u 
n  n  -n  n 

n-1  n=1 

Because  un  is  real  we  must  have  c  =  c  „  as  well  and 
u  n  n 

2 

*  “  iU_t 

u(t)  =  2  Re  T  c  e  u  .  (8.23) 

_  n  n 

n-1 

As  advertised  earlier  our  condition  B*eAtt{/  =  0  for  all  t  e  R+  is 
equivalent  to  the  fourth  component  b(t)  =  0  for  all  t  e  R.  If  we  write 
cR  «  atj,  +  i(5n  and  use  (8.23)  this  means  for  all  t  £  R+  that 


y  (a  sin  y  t  +  0  cos  u  t)  u  w' (L)  =  0  .  (8.24) 

‘‘.n  n  n  n  nn 

n-1 


The  left  hand  side  of  (8.24)  defines  an  almost  periodic  function  and  hence  by 
the  uniqueness  theorem  for  almost  periodic  functions  [19]  the  coefficients 
mush  vanish,  i.e. 

a  p2  w' (L)  =  0  , 

n  n  n 

n  -  1,  2,  . . .  .  (8.25) 

B  u2  w* (L)  =  0  , 

n  n  n 


From  the  spectral  relation  (8.12)  we  know  pn  ^  0.  Furthermore  (8.13)- 
(8.15)  when  coupled  with  w^  *  0  at  x  =  L  force  wn  =  0.  But  this 
contradicts  the  un  being  eigenfunctions  so  we  know  w^(L)  ^  0.  Hence  (8.25) 
inqplies  an  “  “  Or  i.e.  ij>  *  0.  We  can  now  apply  Theorem  7.1  to  conclude 

the  following  result. 


Ih eorem  8.1.  Consider  the  control  system  (8.1)-(8.2).  Then  feedback  control 


yields  a  unique  weak  solution  to  ( 8 . 1 ) ,  (8.2)  for  initial  data  (w^v^a^) 
in  H  -  {(w,v,a,b)  <  H2(0,1)  x  L2(0,1)  xRxR;w  =  w'  =  0  if  x  -  0} .  In 
addition  the  zero  solution  of  (8.1)  is  uniformly  asymptotically  stable  in  H. 

Proof .  Observe  that  u  *  (w,v,a,b),  B*u  =  wfcx(L,t).  So  use  of  Theorem  7.1 
and  the  definition  of  uniform  asymptotic  stability  gives  the  result. 

Notice  that  Theorem  5. 1  also  applies  in  this  case  since  B  is  contact. 

However  as  noted  earlier  Theorem  7.1  yields  the  stronger  result  in  that 
convergence  to  the  zero  equilibrium  solution  is  in  the  strong  rather  than  weak 
topology. 

8.2.  Boundary  control  of  the  NASA  Spacecraft  Laboratory  Control  Experiment  (SCOLE) 
An  exan^sle  where  Theorem  5.1  applies  and  Theorem  7.1  does  not  is  provided 
by  the  stabilization  problem  for  the  NASA  Spacecraft  Laboratory  Control 
Experiment  (SCOLE).  Since  the  details  of  the  mathematical  model  (which  are 
very  lengthy)  are  provided  in  the  papers  of  Balakrishnan  [5]  and  Balakrishnan 
and  Taylor  [6]  we  only  sketch  the  main  ideas  here. 

First  we  quote  from  [5]  for  a  description  of  the  physical  problem. 

"The  physical  apparatus  consists  of  a  softly  supported  antenna 
attached  to  the  space  shuttle  by  a  flexible  beam  like  truss.  Hie  control 
objective  is  to  slew  the  antenna  on  command  within  the  given  accuracy  and 
maintain  stability,  based  on  noisy  sensor  data  and  limited  control 
authority;  allowance  must  be  made  for  random  disturbance.  The  control 
forces  and  torques  are  applied  at  the  shuttle  end  as  well  as  the  antenna 
end  and  in  addition  provision  is  made  for  a  small  number  of  z-axis  proof- 
mass  activators  along  the  beam." 


In  (5]  Balakrishnan  shows  that  in  the  ab sense  of  noise  the  problem  can  be 
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reformulated  in  the  form  (5.6)  where  H  =  (L2(0,L))  x  R  ,  <Au,u>  =  0, 
<Q(u),u>  =  0,  E  =  R1®,  and  B*eAti^  =  0  for  all  t  c  R+  implies  =  0.  So 

the  Corollary  to  Theorem  5. 1  implies  the  zero  solution  is  stable  and 

furthermore  u(t,u0)  +0  as  t  -*■  «  for  all  u0  c  H.  So  the  feedback  control 

law  f(t)  =  K(u(t)}  provides  a  "weakly"  damping  resolution  of  the  SCOLE 
problem  under  the  constraint  f  <  r.  Of  course  since  E  =  R10  and  not  R 

Theorem  7. 1  does  not  apply. 
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